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SECTION  I 


INTRODUCTION 


A.  OBJECTIVE 

Studies  of  soot  formation  and  growth  in  laboratory  combustors  were 
conducted  at  the  Air  Force  Engineering  and  Services  Center  ( AFESC ) ,  Environics 
Division,  using  a  nonintrusive  optical  system  known  as  the  Integrated  Particle 
Sizing  System  (IPSS).  The  behavior  of  this  instrument  and  the  factors 
affecting  its  accuracy  were  not  completely  understood,  causing  the  proper 
operation  of  the  instrument  and  the  interpretation  of  previously  obtained 
results  to  be  questioned.  The  purpose  of  this  research  was  to  investigate 
factors  affecting  instrument  accuracy.  This  insight  allows  the  determination 
of  appropriate  measures  which  optimize  instrument  performance. 

B.  BACKGROUND 

Soot  formation  studies  are  of  interest  because  soot  is  a  product  of 
inefficient  combustion  and  a  good  thermal  radiator.  Soot  deposits  on  the 
engine  liners  of  aircraft  prevent  proper  heat  transfer  and  decrease  the 
lifetime  of  the  liner.  In  addition,  soot  exhaust  emissions  make  the  aircraft 
vulnerable  to  infrared  detector  identification  by  enemy  sensors.  The  presence 
of  soot  can  also  create  a  potential  health  hazard.  Small  soot  particles  (less 
than  1  .0  urn  in  diameter)  can  enter  the  lung  as  they  are  not  easily  deposited  in 
nasal  passages.  It  is  hoped  that  soot  formation  studies  can  result  in 
production  of  a  model  can  which  can  reliably  predict  this  still  uninterpreted 
process  of  combustion. 

Particle  size  distribution  can  be  measured  in  several  different  ways. 
However,  the  latest  trend  is  toward  obtaining  measurements  optically  through 
the  collection  of  radiation  scattered  by  a  particle.  These  optical  diagnostic 
techniques  possess  many  advantages  over  conventional  techniques.  First, 
because  a  focused  beam  of  radiation  serves  as  the  probe,  the  technique 
provides  a  means  of  gathering  information  in  hostile  laboratory  environments 
without  interfering  with  the  flow  field  of  the  combustion  process.  Second, 
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the  spatial  resolution  is  limited  by  the  optics  and  the  wavelength  of  the 
radiation  used,  and  therefore,  for  visible  radiation,  can  be  relatively  small 
(on  the  order  of  microns).  Third,  much  information  can  be  gathered  in  a 
relatively  short  period  of  time  (on  the  order  of  seconds),  thus, leading  to 
good  temporal  resolution  and  statistically  valid  results.  The  IPSS  is  a 
prototype  of  such  an  optical  diagnostic  instrument.  During  preliminary  runs, 
the  IPSS  displayed  unexpected  characteristics  which  are  not  explained  in  the 
system  manual.  This  behavior  prompted  questions  regarding  the  instrument '3 
accuracy  and  initiated  this  study. 

C.  SCOPE 

In  an  attempt  to  answer  questions  about  the  IPSS,  the  theory  governing  the 
instrument's  design  was  studied.  Fhom  this  part  of  the  investigation,  a  list 
of  limitations  imposed  by  theory  was  derived.  The  second  part  of  this 
investigation  studied  the  design  of  the  IPSS.  Work  was  conducted  to  show  that 
the  relationship  between  the  incoming  signal  and  particle  size  does  not  agree 
with  that  found  in  theory.  This  finding  was  verified.  Through  this 
investigation,  a  better  understanding  of  the  IPSS  wa3  gained  and  any 
modifications  to  the  instrument  can  be  done  with  full  understanding  of  their 
impact  on  the  system. 
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SECTION  II 


THE  LIMITATIONS  OF  THE  MIE  THEORY 

The  Mie  Scattering  Theory,  formulated  in  1908,  describes  the 
electromagnetic  field  scattered  by  a  particle  of  arbitrary  size.  The 
expression  for  the  scattered  electromagnetic  field  is  obtained  through  the  use 
of  Maxwell’s  equations  and  the  vector  wave  equation.  In  application,  an 
expression  for  the  irradiance  (incident  energy  per  unit  time  per  unit  area  of 
the  receiver)  is  the  desired  end  result.  In  finding  the  Mie  expression  for 
the  scattered  field,  several  assumptions  were  made  to  simplify  the  problem. 
The  purpose  of  this  section  is  to  point  out  these  assumptions  as  they  appear, 
to  provide  an  awareness  of  the  limitations  imposed  by  the  theory,  and  to  gain 
a  clearer  understanding  of  Mie  theory-based  instrumentation.  Specifically,  a 
foundation  is  formed  which  provides  the  basis  for  the  study  of  factors 
affecting  instrument  accuracy.  In  addition,  the  determination  of  appropriate 
measures  which  optimize  instrumentation  performance  without  violating 
theoretical  considerations  is  also  allowed. 

Geometric  optics  i3  no  longer  valid  when  the  dimensions  of  a  system 
(lenses,  particles,  etc.)  become  comparable  to  or  much  smaller  than  the 
wavelength  of  the  incident  radiation.  Therefore,  for  the  particles  being 
considered,  the  wave  nature  of  light  must  be  examined.  For  this  reason,  a 
brief  review  of  Maxwell's  equations  and  the  boundary  conditions  which  lead  to 
the  Mie  theory  are  included  in  Appendix  A.  Before  the  development  of  the  Mie 
theory,  some  basic  concepts  of  light  interaction  with  matter  and  its 
relationship  to  scattering  must  be  defined. 

A.  SCATTERING  FUNDAMENTALS 

A3  electromagnetic  radiation  travels  through  a  medium,  extinction  occurs; 
that  is,  the  energy  along  the  path  of  the  incident  beam  is  attenuated.  The 
two  processes  responsible  for  extinction  are  scattering  and  absorption.  The 
latter  is  not  the  concern  of  this  paper;  however,  in  simple  terms,  absorption 
is  a  result  of  the  damping  of  the  electromagnetic  wave  by  matter  within  the 
medium.  This,  in  turn,  causes  an  increase  in  the  internal  energy  of  the 
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matter  but  does  not  involve  reradiation  from  the  material.  Scattering,  which 
encompasses  reflection,  refraction,  and  diffraction  occurs  when  the  incident 
radiation  induces  dipole  moments  in  the  particle.  The  dipoles,  oscillating 
with  the  same  frequency  as  the  incident  radiation,  reradiate  the  secondary 
wavelets  in  all  directions.  These  spherical  wavefronts,  which  now  compose  the 
scattered  radiation,  interfere  with  each  other  and  form  an  interference 
pattern  corresponding  to  the  physical  geometry  and  optical  properties  of  the 
particular  scatterer.  Examples  of  such  patterns  are  shown  in  Figures  2  and  3. 

Scattering  can  be  classified  as  elastic  or  inelastic,  dependent  or 
independent,  single  or  multiple.  If  the  scattered  radiation  is  of  the  same 
frequency  as  the  incident  radiation,  the  scattering  is  said  to  be  elastic; 
otherwise  the  scattering  is  Inelastic.  Examples  of  inelastic  scattering  are 
Raman  and  Mandel'stam-Brillouin  scattering.  When  the  presence  of  other 
particles  affect  the  scattering  properties  of  an  individual  particle,  the 
scattering  is  said  to  be  dependent.  Independent  scattering  is  experimentally 
and  theoretically  easier  to  analyze  because  the  total  irradiance  or  intensity 
is  simply  the  sum  of  the  intensities  scattered  by  each  particle.  For 
independent  scattering  to  occur,  the  particles  must  be  "sufficiently  far  from 
each  other"  (Reference  1).  One  general  rule  of  thumb  is  that  the  particles 
must  be  a  distance  of  at  least  three  radii  from  each  other  (Reference  1). 
Reference  2  suggests  that  the  clearances  between  particles  should  be  greater 
than  0.3  wavelength  and  that  the  clearance/diameter  ratio  should  be  greater 
than  0.4.  In  most  cases  the  criteria  for  independent  scattering  is  satisfied 
because  particles  are  usually  separated  by  much  larger  distances  (References  1 
and  2).  The  third  classification,  multiple  scattering,  occurs  when  scattered 
radiation  is  incident  on  another  scattering  surface.  Conversely,  single 
scattering  implies  that  the  scattered  radiation  from  one  particle  propagates 
directly  to  the  detector  without  any  incidence  on  other  particles  that  may  be 
present  in  the  particle  cloud.  Each  scattering  source  scatters  the  radiation 
of  the  original  beam  only.  As  suggested  by  Reference  1  ,  the  prevalence  of 
single  scattering  can  easily  be  determined;  if  doubling  the  concentration  of 
scattering  particles  results  in  twice  the  amount  of  intensity  being  scattered, 
then  single  scattering  is  dominant. 

The  following  pages  show  that  the  scattered  electric  field,  as  derived  by 
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Lorenz  and  later  by  Mie  (References  1,  3.  and  4),  is  a  function  of  the 
magnitude,  wavelength  and  polarization  of  the  incident  field,  the  refractive 
indices  of  both  the  outside  medium  and  the  scattering  medium,  the  size  and 
geometry  of  the  scatterer,  the  distance  from  the  particle,  and  the  direction 
of  scatter. 

B.  DISCUSSION  OF  THE  WAVE  EQUATION 

As  previously  stated,  the  vector  wave  equation  is  the  foundation  of  the 
Mie  theory  and  is  easily  obtained  from  Maxwell's  equations,  provided  certain 
assumptions  are  made  along  the  way.  However,  the  purpose  of  this  discussion 
is  not  to  derive  the  wave  equation,  but  to  show  the  assumptions  leading  to  it, 
making  one  aware  that  other  possibilities  exist  for  a  resulting  fundamental 
equation. 

To  uncouple  the  electric  field  from  the  magnetic  field  in  Maxwell's 
equations,  and  thus  achieve  an  equation  in  terms  of  either  the  electric  field 
or  magnetic  field  only,  the  curl  of  Equations  (A-3)  and  (A-4)  is  taken  and 
Relationships  (A-5),  (A-6),  and  ( A— 7 )  are  substituted  accordingly  to  obtain 

-* 

Vx  V  x  E  =  Vx  (—  ^ 


VxVxH  =  (7xoE)  +  (Vx 


9eE. 

at  1 


(1) 


If  the  electromagnetic  properties,  e,  u,  and  o,  where  e  is  the  electric 
permittivity,  y  is  the  magnetic  permeability,  and  a,  the  conductivity,  are 
assumed  isotropic,  and  if  e  and  u  are  assumed  time  independent,  they  can  be 
placed  outside  the  curl  operations.  Furthermore,  if  the  fields  are  continuous 
in  both  time  and  space,  the  partial  derivative  and  curl  operations  can  be 
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interchanged.  The  result  then  is 


7  x  7  x  E 


■»*&<*  x  H) 


7  x  7  x  H 


o(7  x  E)  +  ej£-(7  x  E) 


Through  another  substitution  of  Equations  (A~3)  through  (A-7)  and  an 
assumption  that  o  is  also  time-independent,  the  equations  become 


2* 

„  _  s  3E  5  E 

7  x  7  x  E  »  -yo^-  ~  UE — 2 

3t 


__  —  ,,  VII 

7  x  7  x  H  -  -yo^-  ' 


With  the  vector  identity 


7  x  (7  x  F)  =  7(7  •  F)  -  7  F, 


the  above  relationships  become 


72I  -  7(7  •  -)  -  yc|£  ♦  ye^-§ 
G  9t2 


72H  -  7(7  •  h  -  yo||  ♦  ye^-J- 
u  3t  3t2 
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Again  using  Maxwell's  equations,  as  well  as  the  assumption  of  isotropic 
electromagnetic  properties,  the  equations  are  written  as 


2-* 
v  E 


Ivp 


uc 


3E 

at 


(6) 


2-» 

V*H 


Assuming  that  both  I  and  H  can  be  written  as  the  time  harmonic  plane  waves, 


E 


t  -i(wt-kz) 
E  e 
o 


H  -  H(e'1<“t~kz>, 


(7) 


where  u  is  the  angular  frequency,  k  is  the  wave  number,  and  z  is  the  direction 
of  propagation,  the  previous  equation  can  be  simplified  to 


V2I  -  -Vp 
e 


( ioiuo 


w2pe)E 


(8) 

V2H  »  (iuiuo  -  ui2ue)H  • 


Now  if  the  gradient  of  the  volumetric  charge  density  (Vp)  is  zero,  and  if  the 
conductivity  is  zero,  or  is  included  in  the  imaginary  part  of  the  refractive 
index  (which  is  related  to  the  relative  permittivity),  the  previous  equations 
become  the  wave  equations: 
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T2!  ♦  «2ueE  -  0 


?■* 
V  H 


2  -*• 
♦  to  peH 


(9) 


One  final  simplification  can  be  made  by  substituting  the  relationship 
u( ue)1/2 -kom,  where  kQ  is  the  wave  number  in  a  vacuum  (-  u)/(uQeo)/2)  and  m  is 
the  refractive  index  (*  (eu/eQuo) ^ )  of  the  medium.  The  wave  equation  can 
then  be  written  in  a  more  familiar  form: 


V2!  ♦  k2m2E  -  0 
o 

2-*  2  2-» 

7  H  +  k  m  H  »  0. 
o 


(10) 


Thus,  the  assumptions  are  as  follows: 

1.  both  e,  u.  and  o  are  isotropic,  time-independent  electromagnetic 
properties, 

2.  E  and  H  can  be  written  as  time-harmonic  plane  waves, 

3.  the  fields  are  spatially  and  temporally  continuous, 

4.  either  the  conductivity,  c,  is  zero  or  is  included  in  the  measurement 
of  the  imaginary  part  of  the  refractive  index. 

The  effect  of  anisotropic  electromagnetic  properties  on  the  scattered 
field  require  that  these  quantities  become  tensors.  The  mathematics  then 
become  much  more  complicated.  These  effects  are  not  investigated  in  this 
report. 

C.  THE  MIE  THEORY 

The  Mie  theory  is  a  mathematically  complicated  expression  which  requires 
the  use  of  a  computer  for  its  calculation.  However,  the  theory  is  a  general 
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one,  unrestricted  by  particle  size.  Other  theories  exist  which  are  less 

cumbersome  to  calculate,  but  their  validities  depend  on  the  ratio  of  particle 

size  to  the  wavelength  of  incident  radiation.  As  a  result,  a  range  of 

particle  size  exists  where  the  Mie  theory  for  particles  of  arbitrary  size 

cannot  be  replaced  by  simpler  theories.  Reference  5  suggests  this  range  to 

be  0.6  A  /im  <  d  <  5  A  /n,  where  n  is  the  real  part  of  the  refractive  index, 
mm 

d  is  the  diameter  of  the  particle,  and  A  ,  the  wavelength  of  the  radiation 

m 

inside  the  particle.  For  diameters  smaller  than  the  lower  bound,  the  Rayleigh 
scattering  theory  is  valid.  For  diameters  larger  than  the  upper  bound,  a 
scattering  theory  for  large  particles  is  appropriate.  Unfortunately,  this 
theory  is  also  popularly  known  as  the  Mie  theory,  although  it  is  a  Mie 
scattering  theory  specifically  for  large  particles.  The  distinction  is  that 
the  Mie  theory  discussed  in  this  report  is  for  particles  of  arbitrary  size. 
Its  formal  name  is  actually  the  Lorenz-Mie  theory;  however,  literature  seldom 
refers  to  it  by  this  name. 

In  addition  to  the  four  assumptions  made  in  developing  the  wave  equation, 
more  assumptions  must  be  made  to  arrive  at  the  Mie  theory: 

1.  the  scattering  is  restricted  to  single,  elastic,  and  independent 
scattering, 

2.  the  particle  is  spherical, 

3.  the  particle  material  is  homogeneous, 

M.  the  magnetic  permeability  of  the  particle  is  equal  to  the  magnetic 
permeability  of  the  outside  medium, 

5.  the  relative  refractive  index — that  is,  the  ratio  of  the  refractive 
index  of  the  particle  to  that  of  the  outside  medium — mus-t  be  known, 

6.  no  free  surface  charge  is  present  on  the  particle, 

7.  the  incident  radiation  is  plane  wave,  and 

8.  the  scattered  radiation  is  measured  far  from  the  point  of  scatter  so 
that  the  longitudinal  electric  field  is  negligible. 

The  schematic  of  the  problem  is  shown  in  Figure  1.  The  azimuthal 

angle  Q  is  defined  as  the  angle  formed  by  the  propagation  vectors  k.  (the 

direction  of  propagation  of  the  incident  radiation)  and  k  (the  direction  of 

s 

propagation  of  the  scattered  radiation).  The  plane  containing  this  angle  is 
called  the  plane  of  scatter.  Once  this  plane  is  defined,  the  incident 
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electric  field  can  be  decomposed  into  two  components:  Ej^  and  Ej  which  are 

perpendicular  and  parallel  to  the  plane,  respectively.  Note  that  the 

magnitudes  of  these  components  are  a  function  of  the  polar  angle  <(»  which  is 

formed  by  the  x-axis  and  the  projection  of  k  on  the  xy-plane.  Further  note 

s 

that  since  $  depends  on  the  orientation  of  the  plane  of  scatter,  the 
magnitudes  of  the  perpendicular  and  parallel  components  of  the  incident 
radiation  also  vary  with  the  orientation  of  the  plane  of  scatter. 

The  desired  result  of  the  Mie  theory  is  an  expression  for  the  electric 
field  scattered  by  the  particle.  In  its  simplest  and  most  general  form,  the 
equation  can  be  written  as 

(E|.s\  gik(r-z) 

El,s/  '  “ikr 

where  [S]  is  the  amplitude  scattering  matrix  consisting  of  complex  numbers. 
As  will  be  seen,  this  matrix  is  governed  by  the  properties  and  geometry  of  the 
particle,  and  the  azimuthal  angle  of  scatter. 


S„(0)  So(0) 


s4(9)  s1(e)|\Ei  ±(d>) 


i  KM 

J  Vi  i*7 


(n) 


In  application,  the  intensity,  I,  of  the  electric  field  i3  measured. 
Thus,  the  incident  and  scattered  electromagnetic  radiation  are  frequently 
described  by  the  Stokes  Parameters  I,  Q,  U,  and  V,  which  are  defined  below: 


1  ’  S3  CE|E|*  *  EiEi*] 

Q  '  S3  CE|E|*  -  EiEi*] 

(12) 

"  '  S3  tE|EJ*  *  EJ_E|*1 
V  ■  755  CE|Ei*  -  ElT]- 
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With  the  aid  of  Equation  (11)  and  following  the  notation  of  Reference  4,  the 
scattered  radiation  can  be  characterized  in  terms  of  the  incident  radiation  by 


where 

S11  *  I  (ISJ2  +  lS2!2  +  lS3|2  *  lSJ2) 

si2  =  i  (IS2I2  ~  lsJ2  *  ls4l2  -  Is3I2) 

513  “  Re{s2S3*  +  SlV} 

514  -  Im{S2S3*  -  S^*} 

521  *  \  (|S2|2  '  lSl!2  ‘  I S4 I ^  +  Is3!2) 

522  “  \  (lS2i2  +  lSJ2  "  I S4 1 2  "  iS3f2) 

523  "  RetS2S3*  -  W) 

524  -  Im{S2S3»  +  S^*} 

S31  -  Re(S2S4»  ♦  S1S3»} 


(13) 
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S33  -  Re[S,S2*  *  S3S4*} 
S3J<  -  ImlSjjS^  *  S4S3«} 
S„1  -  Im(S2*S4  ♦  S3»S1} 

542  -  Im{S2«S4  -  S^S,} 

543  -  Im{S^S2*  -  S3S4»} 
s44  -  Re{S1S2*  -  S3S4»}  . 


(H) 


Therefore,  the  scattered  electric  field  and  scattered  intensity  can  be 
described  if  S^e),  S2(e),  S3(0),  and  S^C 0)  are  found.  These  values  are 
found  by  solving  the  wave  equation.  However,  the  search  i3  simplified  by 
knowing  that  for  homogeneous  spherical  particles  S3(e)  -  S4(0)  »  0.  Thu3,  the 
scattered  electric  field  is 


E 

s 


ik(r-z) 


♦  S2(  9)E|  ^ ( <^ ) a j  ^  , 


(15) 


and  the  scattered  intensity  is  simply 


[IVe)l2 «!.,  ‘ 


(16) 


To  begin  the  search,  the  wave  equation  must  be  written  in  spherical 
coordinates  since  spherical  particles  are  considered.  Unlike  the  position 

A  •* 

vectors  in  the  Cartesian  coordinate  system,  the  position  vectors  a^,  ag, 

and  a.  are  not  fixed.  Thus,  the  components  of  E  separately  will  not  satisfy 
<P 

the  scalar  wave  equation 
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T2  ip  +  k^4>— 0 . 


(17) 


However,  as  will  be  seen,  the  solution  of  the  scalar  wave  equation  is  useful 
in  finding  a  solution  to  the  vector  wave  equation.  The  solutions  of  ip  can  be 
3hown  to  be 

o 

ip  .  -  cosi$  P  (cos0)Z  (kr) 
efcn  n  n 

^  ip  .  -  sin£$  P*'(cos0)Z  (kr)  . 
o*.n  n  n 


where  subscripts  e  and  o  refer  to  the  even  and  odd  solutions  of  the 
function  ip, 
integers  n  i  l  i  0, 


Pn  is  the  associated  Legendre  polynomial  of  the  first  kind  of  degree 


n  and  order  1, 


is  one  of  the  spherical  Bessel  functions  listed  below: 
Vo 


Jn(kr)  - 

(2krJ  Jn-*V2(kr) 

spherical 

kind. 

-  v2 

yn(kr)  = 

!■ 

1 

H< 

D 

X 

spherical 

kind, 

h(1)  (kr) 
n 

*  jn(kr)  *  iyno<r) 

spherical 

kind, 

h^2)  (kr) 

-  Jfi(kr)  ■  iyn(kr) 

spherical 

kind. 

The  transformation  which  permits  the  solution  of  the  vector  wave  equation 
via  the  scalar  wave  equation  is 


M  »  7  x  (rip) . 


09) 


where  r  »  ra^.  The  scalar  wave  equation,  showns  that  this  transformation 
satisfies  the  vector  wave  equation  because 


2-* 

7*M  ♦ 


k2M 


**  x  [r(72<p 


+  k2ip)  ] 


0. 


(20) 
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Also,  like  the  electric  and  magnetic  field,  the  divergence  of  the  vector  M  is 
zero  because  the  divergence  of  a  curl  of  any  vector  is  always  zero. 
Furthermore,  M  is  coupled  to  another  vector  i$  by 

V  x  N  -  kM 

(21) 

V  x  M  -  kN  . 

-♦ 

The  two  solutions  of  M  can  be  expressed  in  terms  of  the  even  and  odd 
functions  of  41: 


M  „  ■  7  x  n|»  . 

edn  redn 


M  „  -  7  x  ri|»  .  . 

odn  odn 


(22) 


-+ 

of  N  can  be  expressed  in  terms  of  and 


(23) 


Reference  4  defines  M  and  N  as  vector  harmonics  and  point 'out  that  both 
have  the  qualities  of  an  electromagnetic  field — that  is,  both  satisfy  the 
vector  wave  equation,  both  have  a  zero  divergence,  and  both  have  curls  which 
are  proportional  to  each  other.  Since  these  vector  harmonics  are  mutually 
orthogonal,  they  are  useful  in  obtaining  an  expression  for  the  incident 
electric  field.  The  electric  field  can  now  be  written  as  an  expansion  in 
terms  of  the  vector  harmonics: 


Likewise,  the  two  solutions 


w 


7  x  M 


edn 


edn  k 


7  x  M 


odn 


odn  k 


.1  LAB  „  M  + 
d-o  n-d  edn  edn 


B  „  M  ♦ 
odn  odn 


A  N  * 
edn  edn 


AodnNodn^ ’ 


where 


(24) 
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M 


eln 


pj(cose)  „  dP£(cO30) 

-lsini*  sine  '  Zn(kr)ae  -  C03**  ^de - zn(kr)a4 


P^(cos0)  „  P^"  ( GOS0) 

Kin  *  lCOSl*  — n—  zn(kr)a0  -  sin**  - di -  Zn(kr)ad 


„  l  Zn(kp)  ' 

Nein  "  n(n+1)c°si^  Pn(cos0)  — -■■■-  ap 


+  cost# 


dP*'(cos0)  ,  d[krZ  (kr)]  . 
n  i  n 


d0 


kr  d(kr )  a0 


P  (co30)  .  d[krZ  (kr)3  ^ 
nHe  kr  dTkrl  % 


„  l  zn(kr)  - 

N  „  »  n(n+1  )3lnl<>  P  (cos0)  — : -  a 

oUn  n  kr  r 


(25) 


P*"( cos0)  d[krZ  (kr)] 

+  sinH  -Sjg -  kf  “~dTk"rV  3 


0 


P  (cos0)  d[krZ  (kr)]  « 

+  2,cos2.<J)  -  ■— —  - -tt; — r — 

sin0  kr  d(kr)  <j 


A  „  ,  A  0  ,  B  .  ,  and  B  „  are  coefficients  whose  values  are  derived  from 
eln  oin  ein  oln 

orthogonality  conditions  of  the  inner  products  of  (E^,  Mgin)»  ^Mej_n»  Mej,'n'^’ 

<B1>  W-  <W  <Er  W'  (IW  W’-  <Er  an“ 

(Notn*  Nol'n,)' 

For  the  incident  linearly  polarized  plane  wave, 


e,  .  e 

i  o  x 


E  i(krco30  uit) (ain0cos$a  +  cos0cos$a.  -  sin<|>a,), 
o  r  O  q> 


(26) 
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J. 


the  resulting  values  of  these  coefficients  are  as  follows: 


A  *  A 
Aetn  eln 


-iE 


o 


in(2n+1 ) 
n(n+1 ) 


Bo1n’ 


Aoln 


0, 


(27) 


B 


eln 


0, 


B  „  -  B  , 

oln  oln 


ln(2n+1) 
o  n(n+1) 


Only  nonzero  coefficients  are  obtained  from  orthogonality  conditions 
when  4*1.  Thus,  the  expansion  of  the  incident  electric  field  in  spherical 
harmonics  is 


E  Z.  i 
o  n-1 


2n+1 
n(n+1 ) 


CM 


oln 


iN  ,  ] 
eln 


(28) 


Since  the  magnetic  field  is  coupled  to  the  electric  field  by  Faraday's 
law,  H  -  -  i(V  x  E)/uu,  the  incident  magnetic  field  can  be  found  by  applying 
Equations  (21)  and  (23)  to  the  expression  for  the  incident  electric  field  to 
obtain 


k_  E  £  tn  2n+1 
up  o  n-1  n(n+1 ) 


^oln5 


(29) 


Expressions  for  the  internal  and  scattered  field  follow  from  these 
expansions.  Following  the  notation  of  Reference  4,  the  field  within  the 
particle  is  written  as 


„  2  .n  2n+1  ,2 

E  Z,  i  - rr  (c  M  ,  -  id  N  ,  ) 

o  n-1  n(n+1 )  noin  n  el n 
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(30) 


-  !<£_  v  ?  in  2n'f1  (a  5 
wu  o  n-1  n(n+1)  n  eln 


ic  N  ,  ). 
n  oln 


Likewise,  the  scattered  field  is  expressed  by 


K  ■  Eo  nil  *"  sfSiTT  ‘-Vo.n.n  *  lV„n,h> 


k  « 

—  E  E,  i 
o)ii  o  n-1 


n  2n+1  ,  j» 

n(n+1 )  anMe1  n,h  + 


ib„NMn  J. 
n  o i n, h 


(3D 


The  coefficients  an,  bn,  cn,  dn  are  found  by  applying  boundary  conditions  for 
a  perfect  dielectric  (Appendix  A  8a-b)  to  the  fields  tangential  to  the 
particle's  surface,  namely  (E_,  H.)  and  (Ex,  H  ).  Thus,  it  is  assumed  that 

9  u  cp 

no  net  free  charge  is  present  on  the  particle's  surface.  If  the  magnetic 
permeability  of  the  particle  is  also  assumed  to  be  equivalent  to  that  of  the 
surrounding  medium,  the  expression  for  the  coefficients  can  be  simplified  to 


'*n(y)’*n(x)  “  mVyH;(x) 

3n  ”  ^y^n(x)  ~  m\pn(y)^(x) 

m^(y)^n(x)  -  ♦  (y)^(x) 
bn  =  m^ ( y ) £n ( x )  -  <l/n(y)Z^(x) 

mi jin(x)E^(x)  -  m5n(x)i|^(x) 

°n  ”  *n(y)^(x)  -  m£n ( x ) ^ ( y ) 

m^n(x)^(x)  -  mCn(x)iii^(x) 
dn  ”  mil»n(y)^(x)  -  £n(x)t|^(y)  ’ 


where  i|»n(p)  -  pJn (  p ) 

^n(p)  -  ^(p) 

x  -  k.a  -  2ira/X  , 
°  o 

y  -  ka, 


18 


k  =  mk0, 

a  «  radius  of  the  particle. 

The  vector  harmonics  are  written  as 

Ln  -  -sin*  n  (cose)Z  (k  r)a  -  cos*  x  (cos0)Z  (k  r)a, 
eiri  n  n  ra  w  n  n  ni  ^ 


oln 


cos*  II  (cos9)Z  (k  r)a.  -  sin* 
n  n  m  9  r 


x  ( cos 0 ) Z  (k  r)a, 
n  n  ro  d 


„  Z  (kr)  . 

N  .  =  n(n+1)cos*  sin0  n  (cos@)  -n  —  — —  a 

ein  n  k  r  r 


.  d[k  r  Z  (kr)]  „ 

+  COS*  X  (COS0)  - - J7r - r -  a 

n  kr  d(k  r)  0 

m 


1  dtkmr  Zn(kmr)^  - 

-  sin*  n  (cose)  ,n-  ” -  a 

n  k  r  d(k  r)  * 

mm 

-  Zn(kmr)  - 

NQln  =  n(n+l )sin*  sine  Iln(cose)  ■  ■  k  ” —  a^ 


.  d[k  r  Z  (k  r) ]  » 

+  sin*  x  (cose)  - -  - 2-r- — 2-r — -  a 

n  k  r  d(k  r)  9 

m  m 


,  d[k  r  Z  (k  r)]  . 

+  COS*  n  (003  9)  r— - -T, ,  -  rS -  a 

n  k  r  d(k  r)  d 

m  m 


where  n  (cose) 
n 


dP  (cose) 
n 


p'  ( cose) 
n 


sine  d(cos9)  sine 


x  (cose) 
n 


.  dP  (cos9) 
a  n _ 

de  d(cos9) 


de  p;,(co39)' 


km  *  k  of  the  particular  medium  where  the  electromagnetic 
propagation  is  considered, 
r  »  distance  from  the  particle  center. 


(33) 
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For  considering  the  electromagnetic  fields  incident  on  and  within  the 

particle,  the  appropriate  spherical  Bessel  function  used  in  the  expressions 

for  M  and  N  is  J^^r)  because  the  fields  are  finite  at  the  origin  (the  center 

of  the  particle).  For  the  scattered  electromagnetic  field  the  appropriate 

spherical  Bessel  function  describing  an  outgoing  spherical  wavefront  is  the 

spherical  Hankel  function  of  the  first  kind,  h^(k  r):  thus,  the  subscript  h 

n  m 

is  present  in  the  expression  for  the  scattered  field. 


Furthermore,  for  large  distances  from  the  particle  (kmr  >>  n2),  the 
spherical  Hankel  function  is  approximated  as 


(kr) 


,  ..n  ikr 
(-i)  e 


(34) 


The  use  of  this  approximation  leads  to  a  (kr)2  term  in  the  denominator  of  the 
radial  component  of  the  scattered  electromagnetic  wave.  Therefore,  at  large 
distances  from  the  particle,  the  scattered  radial  field  can  then  be  neglected. 


By  substituting  all  the  appropriate  functions  into  the  expression  for  the 
scattered  electric  field,  simplified  expressions  for  the  components  are 
obtained: 

Ik  r 

n  .  m 

r  f  -iwt  oC  e  «  (2n+1 )  ,  . 

Es,e  =  9,|  “  e  -ik  r  n»1  n(n+1)  anTn  +  bnnn 

”  (35) 

ik  r 

p  ,  p  =  -i»t  Eo3in*  6  m  -  (Zn+i)  .  _  +  .  . 

Es,i>  s,|  -ik  r  n-1  n(n+l)  an^n  bnTn 

fTJ 


However,  early  in  the  derivation,  the  polarization  of  the  incident  radiation 
wa3  assumed  to  be  along  the  x-axis.  As  a  result,  the  parallel  and 
perpendicular  components  of  the  incident  electric  field  can  be  written  as 


Eocos<j>  e 


i(kz-uit) 


EQsin$  e 


i(kz-wt) 


(36) 
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With  these  equations,  the  final  expression  of  the  scattered  electric  field  is 
obtained: 
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Thus,  for  spherical  particles,  the  search  for  the  elements  of  the  matrix 
[S]  is  complete.  By  comparing  the  final  result  with  Equation  (11),  and  by 
remembering  S^(0)  -  3^(0)  -  0  for  spherical  particles,  it  can  be  seen  that 


s1  (0) 


«  2n+1  .  _  .  . 

Z,  —> — rr  (an  +  b  t  ) 
n-1  n(n+1 )  n  n  n  n 


(38) 
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Consider,  for  example,  the  case  of  scattering  in  the  4>  -  90°  plane  (yz- 
plane).  In  this  circumstance,  the  incident  radiation  is  perpendicularly 
polarized  with  respect  to  the  plane  of  scatter.  From  the  previous  equation, 
the  scattered  intensity  is  then 


I 


s 


.2  2 
k  r 


(39) 


In  Figures  2  and  3.  the  common  logarithm  of  this  scattered  intensity 

normalized  by  the  factor,  IQ/k^r^,  is  plotted  against  the  azimuthal  angle  0 

for  various  particle  sizes  and  two  refractive  indices.  In  these  diagrams, 

light  is  incident  at  0  »  180°;  thu3,  backscatter  occurs  for  0  -  180°  and 

0 

forward  scatter  occurs  for  8  »  0  .  Figure  2  is  obtained  by  assuming  a  real 
refractive  index  while  Figure  3  shows  the  intensity  profile  when  the 
refractive  index  is  a  complex  number.  As  can  be  seen  from  comparing  Figures 
2c  and  3c,  the  presence  of  the  imaginary  part  in  the  refractive  index  is 
responsible  for  the  damping  or  attenuation  of  the  electromagnetic  wave. 
Through  a  comparison  of  the  intensity  profiles  for  each  particle  it  can  also 
be  seen  that  as  the  particle  size  increases,  the  forward  scatter  lobe  becomes 
narrower  and  the  magnitude  of  the  forward  scattered  light  increases. 
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In  some  circumstances  (i.e.,  particle  sizing  measurements  in  combustion 
environments),  the  incident  intensity  may  be  unknown  and  difficult  to 
measure.  This  unknown  parameter  presents  a  problem  if  other  quantities  such 
as  particle  size  are  also  unknown.  However,  this  problem  can  be  circumvented 
by  taking  the  ratio  of  the  scattered  intensities  collected  at  two  different 
angles  of  scatter.  Mathematically,  for  linearly  polarized  light  perpendicular 
to  the  plane  of  scatter  with  $  fixed  at  90°,  the  result  is  an  equation 
relating  the  intensity  ratio  to  particle  size: 


X  , 

ratio  's.B 


(40) 


Plots  of  intensity  ratio  versus  particle  size  for  a  fixed  angle  pair  and 
refractive  index  are  shown  in  Figures  4a  and  4b.  Note  that  for  a  complex 
refractive  index,  the  curve  is  relatively  well-behaved;  monotonically 
decreasing  for  increasing  particle  size.  However,  this  behavior  is  not  seen 
for  a  real  refractive  index.  Although  a  given  particle  size  corresponds  to  a 
specific  intensity  ratio,  the  converse  is  not  true.  This  multiplicity  does 
not  pose  a  serious  problem  in  determining  particle  size  if  it  is  remembered 
that  at  near  forward  scattering  angles,  the  absolute  scattered  intensity 
increases  with  increasing  particle  size.  As  a  result,  provisions  must  be  made 
in  the  electronics  of  Mie  theory-based  instrumentation  to  allow  the 
discrimination  of  different  size  particles  having  the  same  intensity  ratio. 
One  specific  size  discrimination  scheme  is  discussed  in  the  next  section. 

Many  computer  codes  (References  4,  6,  and  7)  are  available  which  calculate 
the  scattered  intensity.  However,  it  is  important  to  note  that  the 
expressions  shown  are  only  point  expressions  for  the  scattered  field.  In 
practical  situations,  the  scattered  radiation  is  collected  over  a  finite  solid 
angle.  Mathematically,  this  measurement  translates  into  the  problem  of 
finding  the  total  scattered  intensity  at  various  planes  of  scatter.  In  other 
words,  the  point  values  must  be  integrated  over  the  desired  solid  angle.  As 
can  be  seen  by  the  complexity  of  the  Mie  expression,  a  numerical  integration 
scheme  must  be  implemented. 
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SECTION  III 


THE  INTEGRATED  PARTICLE  SIZING  SYSTEM 

The  instrumentation  used  for  nonintrusive  measurement  of  particle 
diameters  in  hostile  environments  is  composed  of  an  optical  system,  a  data 
processing  unit,  and  a  data  management  system.  To  test  and  calibrate  the 
system,  a  particle  generation  system  for  known  particle  sizes  is  also 
necessary.  In  the  following  pages,  each  component  is  discussed,  as  well  as 
the  calibration  and  practical  application  of  the  instrument. 

A.  OPTICAL  SYSTEM 


The  Integrated  Particle  Sizing  System  (SDL  Model  IPSS-200),  designed  by 
Spectron  Development  Laborator ies,  Inc.  (Costa  Mesa,  California),  to  measure 
the  formation,  growth,  and  destruction  of  soot  particulates  in  flames,  has  a 
large  measurement  range  with  the  incorporation  of  two  techniques  for  particle 
sizing:  scattered  intensity  ratioing  for  particle  diameters  between  0.08  and 
3.0  um  and  particle  sizing  interferometry  for  particle  diameters  between  2.5 
and  45  um.  Particle-sizing  interferometry  is  an  extension  of  the  laser  Doppler 
velocimeter  technique.  Scattered  intensity  ratioing,  which  forms  the  ratio  of 
the  scattered  light  collected  at  two  different  angles,  can  be  further 
subdivided  into  large-angle  and  small-angle  ratioing  techniques.  The  main 
difference  between  these  two  techniques  is  the  configuration  of  the  receiving 
optics.  As  described  by  Spectron,  the  small-angle  technique  (8°/2.5°,  10°/5°, 
and  12°/5°)  uses  annular  masks  to  define  the  separate  angles.  Large- 
angle  ratioing  (40°/20°,  50°/20°,  and  60°/20°)  uses  a  separate  set  of  lenses 
for  each  angle  of  the  angle  pair.  The  dynamic  range,  as  reported  by  Spectron, 
for  the  complex  refractive  index  of  soot,  1 . 57— i 0.56,  and  the  six  angle  pairs 
used  for  the  small-angle  and  large-angle  ratioing  technique  is  given  in  Table 
1 .  Note  that  for  the  study  of  soot  formation  from  gas  turbine  engines 
(typically,  the  soot  particles  average  0.15  um  in  diameter),  the  large-angle 
technique  is  of  particular  interest.  Thus,  the  large-angle  configuration  will 
be  described. 

A  schematic  of  the  IPSS  is  shown  in  Figure  5.  The  incident  radiation  is 
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TABLE  1.  DYNAMIC  RANGE  FOR  SOOT  AS  PROPOSED  BY  SPECTRON. 
For  1.57  -  10.56  (soot) 


ANGLE  PAIR (degrees) 

MINIMUM  DIAMETEROm) 

MAXIMUM  DIAMETEROm) 

8°/2.5° 

0.45 

2.95 

1 0°/5° 

0.40 

2.44 

12°/5° 

0.30 

1.97 

40°/20® 

0.12 

0.57 

50°/20° 

0.09 

0.45 

60°/20° 

0.08 

0.37 
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Figure  5.  Schematic  of  the  IPSS 


supplied  by  the  0.488  ym  line  of  a  2-Watt  Lexel  Argon  Ion  Laser  (Model  95) 
operating  in  the  multiline  mode.  The  various  lines  are  separated  by  a 
dispersion  prism;  however,  only  the  0.488  um  line  is  allowed  to  pass  through 
the  system.  After  reflecting  from  mirrors  Ml,  M2,  M3,  and  M4,  the  beam  is 
expanded  by  a  variable  beam  expander  comprised  by  lenses  LI  and  L2,  which, 
according  to  Spectron,  has  the  potential  of  increasing  the  beam  diameter  over 
a  range  of  1.7  to  8.3  times  its  initial  size.  After  expansion,  the  beam  is 
directed  toward  the  dichroic  mirror  which  ensures  that  the  electric  field  is 
linearly  polarized.  The  beam  is  then  focused  by  an  800  mm  lens  (L3).  The 
tightly  focused  beam  provides  a  small  sampling  volume,  good  spatial 
resolution,  a  high  signal-to-noise  ratio,  and  a  method  of  counting  particles 
if  only  one  particle  at  a  time  is  in  the  probe  volume  for  each  measurement. 
Thus,  if  the  waist  diameter  of  the  focused  beam  is  theoretically  related  to 
the  beam  diameter  in  the  equation, 

dua)o.  -  4Xf/ird,  (41 ) 

waist  laser 

where  dwai3t  -  diameter  of  the  beam  at  the  waist, 
dlaser  ”  diameter  of  beam  before  focusing, 
f  »  focal  length  of  focusing  lens, 

X  *  wavelength  of  the  radiation, 

the  purpose  of  beam  expansion  before  focusing  is  realized. 

The  radiation  scattered  by  particles  entering  the  probe  volume  is 
collected  at  one  of  the  previously  determined  angle  pairs  by  lenses  L4,  L5, 
L6,  and  L7  and  focused  onto  photomultiplier  tubes  PMT  A  and  PMT  B  (RCA  Model 
8575).  Included  in  the  PMT  housings  are  bandpass  filters  which  have  a 
bandwidth  of  approximately  2  nm  and  a  transmittance  of  approximately  20 
percent  (Figure  6).  Power  is  supplied  to  the  PMTs  by  a  triple-high-voltage 
power  supply  (SDL  HV-1000)  which  can  operate  three  PMTs  and  is  capable  of 
generating  voltages  up  to  -1000  volts  dc  at  4  mA.  As  with  all  photomultiplier 
tubes,  the  amount  of  light  observed  at  the  aperture  of  the  tube  is  directly 
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Transmittance  (100%  Full  Scale) 


Wavelength  (um) 


Figure  6.  Percent  Transmittance  versus 
Wavelength  Measurements  for 
PMT  Filters. 


30 


proportional  to  the  negative  output  current  when  operating  in  the  tube's 
linear  current  range.  The  resulting  signals  from  the  two  PMTs  are  then  sent 
to  the  processing  units. 

B.  DATA  PROCESSING  SYSTEM 

The  PMT  output  currents  are  first  passed  through  a  logarithmic  amplifier 
(SDL  LA-1000)  where  they  are  converted  to  positive  voltages  and  are  scaled  by 
a  factor  of  2  volts  per  decade  current  with  an  input  current  of  -1  mA 
producing  a  10  volt  peak  output  pulse.  In  equation  form,  the  relationship  can 
be  written  as 

vs  *  21og(i3/i0) 

where  vg  =*  output  voltage  of  the  signal,  (42) 

i3  -  signal  current  in  amperes, 


The  voltage  signals  are  then  directed  to  the  ratio  processor  (SDL  RP-1001). 
Within  the  processor,  several  operations  take  place.  First,  the  signals  are 
passed  through  a  peak  detection  circuit.  A  peak  is  said  to  occur  when  the 
large-angle  signal  goes  above  the  1 .5~volt  threshold  value  and  then  decays  to 
a  value  of  1  volt  below  its  corresponding  maximum  voltage.  The  detection 
circuit  is  reset  when  the  large-angle  signal  goes  below  the  1-volt  threshold 
level.  The  dc  levels  of  the  two  peaks  are  then  sent  to  an  analog  sub  tractor 
circuit  which  subtracts  the  large-angle  signal  (A)  from  the  small-angle  signal 
(B).  The  resulting  signal  becomes  the  common  logarithm  of  the  intensity 
ratio,  since 

VB  -  VA  =  -21og(iA/iB).  (43) 

The  analog  subtractor  circuit  has  an  output  ranging  from  0  volt  to  2  volts  for 
the  ratio  iA/iB  of  1  to  0.1,  respectively.  The  range  of  the  subtractor  output 
is  increased  to  0  to  10  volt3  when  the  signal  is  passed  through  an  amplifier 
with  a  gain  of  5.  Upon  amplification,  the  output  travels  to  an  analog/digital 
converter  where  it  is  converted  to  an  eight-bit  binary  number.  The  validity 
of  the  data  word  i3  internally  checked  in  the  processor  by  examining  the 
eighth  bit.  Thus,  the  eighth  bit  is  used  as  a  status  bit.  As  outlined  by 
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4-volt  Threshold 


1.5  V  \  (c). 


time 


time 


Figure  7.  Examples  of  the  Restrictions  of  Valid  and  Invalid  Data: 

(a)  A  possible  signal,  but  signal  does  not  deteriorate 
to  1  volt  below  the  maximum  before  the  next  signal 
arrives.  This  signal  is  not  recoqnized  by  the 
processor . 

(b)  A  valid  data  point. 

(c)  The  peak  signal  is  less  than  the  large-angle  signal 
of  4  volts  (invalid  data) . 

(d)  The  pulse  width  is  less  than  2  ysec  (invalid  data) . 

(e)  The  ratio  of  the  large-angle  signal  to  the  small- 
angle  signal  is  greater  than  1.0  (invalid  data). 

(f)  The  peak  of  the  small-angle  signal  is  greater  than 
the  small-angle  threshold  value  of  10  volts  (invalid 
data 


Spectron,  a  data  word  is  considered  valid  when  the  following  three  criteria 
are  met: 


1.  The  input  signal  level  from  the  small  angle  must  be  less  than  1C 
volts.  In  terms  of  current,  the  output  of  PMT  B  must  be  less  than  1 
mA. 

2.  The  input  signal  level  from  the  large  angle  must  be  greater  than  a 
minimum  of  +4  volts  for  high  speed  operation.  In  terms  of  current, 
the  output  of  PMT  A  must  be  greater  than  1  uA. 

3.  The  width  of  the  input  signal  pulse  from  the  large  angle  must  be 
greater  than  a  minimum  of  2  usee  for  high  speed  operation. 


Figure  7  demonstrates  these  restrictions.  As  a  consequence  of  Restrictions  1 
and  2,  a  particle  size  outside  the  range  of  a  specified  angle  pair  should  not 
be  considered  as  valid  data  although  that  3ize  may  have  an  intensity  ratio 
equal  to  that  of  a  particle  inside  the  angle  pair's  range.  For  example, 
although  latex  particles  of  0.456  ym  and  1.39  pm  have  the  same  intensity  ratio 
at  an  angle  pair  of  40°/20°  (Figure  8),  the  absolute  scattered  intensity  for 
the  larger  size  particle  is  about  10  times  that  of  the  smaller  particle 
(Figure  9).  Therefore,  if  the  electronics  of  the  system  are  properly 
adjusted,  the  signal  from  the  large  particle  is  rejected. 


C.  DATA  MANAGEMENT  SYSTEM 


The  digitized  data  is  transferred  to  an  Apple  II  microcomputer  which,  with 
the  help  of  the  software  package  (SDL  HG-3001).,  stores  the  real-time  data  in 
memory  and  displays  a  histogram.  These  histogram  results  are  also  stored  in 
tabular  form.  An  example  of  the  data  is  shown  in  Figure  10.  The  intensity 
ratio  range  of  0.1  to  1.0  is  divided  into  64  subintervals  or  bins  according  to 
the  equation 


where  N 


T  in-(N  +  BlN#)/64 

^bin  ” 

!0,  for  maximum  bin  ratio 
1/2,  for  bin  center  ratio 
1,  for  minimum  bin  ratio  . 


(44) 
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These  ranges  are  shown  in  Table  2.  Particle  sizes  are  correlated  with  each 
bin  intensity  ratio  range  by  use  of  bin  tables  provided  in  the  IPSS  software 
package.  Note  that  only  62  bins  are  shown  in  the  histogram.  Bins  1  through 
62  are  treated  as  valid  data  and  are  displayed.  However,  since  the  ratio 
processor  will  accept  only  those  intensity  ratios  which  are  less  than  1.0  but 
greater  than  0.1,  Bin  0  which  contains  any  ratio  greater  than  or  equal  to  1.0, 
and  Bin  63,  which  contains  any  ratio  less  than  or  equal  to  0.1,  are  treated  as 
invalid  data.  The  display  is  limited  to  152  data  points  (counts)  in  each 
bin.  Data  acquisition  is  terminated  when  either  a  bin  becomes  full,  the  total 
number  of  data  samples  stored  in  the  histogram  exceeds  a  value  set  by  the 
operator  (limited  by  the  values  101  to  10^),  or  the  operator  terminates  the 
procedure.  Stored  within  the  table  will  be  the  maximum  raw  count  (the  maximum 
count  in  one  bin),  the  total  raw  count  (the  sum  of  all  the  counts  in  the  valid 
bins  1  through  62),  and  the  sampling  time.  If  desired,  the  linear,  surface, 
volume,  and  Sauter  mean  (also  known  as  the  mean  volume-surface)  diameters,  as 
well  as  the  standard  deviation,  can  be  calculated.  The  sample  output  shown  in 
Figure  10  displays  the  bin  center  intensity  ratios,  its  corresponding  bin 
center  diameter,  and  the  number  of  counts  in  each  bin. 

Observation  of  a  Tektronix  dual-trace  oscilloscope  (Model  434)  and  the 
ratio  display  and  data  rate  display  on  the  ratio  processor  panel  aid  in  the 
real-time  monitoring  of  these  signals,  as  well  as  in  determining  whether  or 
not  the  above  criteria  are  met.  The  oscilloscope  monitors  the  logarithmic 
amplifier  output  and  displays  the  voltage  signal  of  each  angle.  The  intensity 
ratio  display  shows  the  previously  calculated  ratio  at  a  1  Hz  rate.  The 
presence  of  large  or  small  particles  out  of  the  range  of  detection  is  signi¬ 
fied  by  a  ratio  display  of  0.1  or  0.99,  respectively.  The  data  rate  display, 
also  updated  at  a  1  Hz  rate,  indicates  in  kHz  the  number  of  signals  read  in 
the  previous  time  interval.  A  low  data  rate  can  indicate  that  either  a  low- 
or  high-number  density  is  present.  In  the  latter  case,  the  large-angle  signal 
does  not  deteriorate  to  the  1-volt  threshold  level  before  another  signal  peak 
occurs.  By  monitoring  the  signal  on  the  oscilloscope,  the  operator  can  deter¬ 
mine  if  either  situation  is  prevalent — square  pulses  in  conjunction  with  a  low 
data  rate  suggest  that  a  large  number  of  particles  are  present  in  the  probe 
volume.  A  maximum  theoretical  data  rate  of  100  kHz  has  been  approximated  by 
Spectron  by  using  Restrictions  1  and  3  outlined  previously. 
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TABLE  2.  LIST  OF  INTENSITY  RATIO  BIN  MAXIMUMS, 
MINIMU MS ,  AND  CENTERS  FOR  EACH  BIN. 


BIN  RATIO  RANGES 


BIN 


BIN 

WIDTH 

MAXIMUM 

MINIMUM 

CENTER 

33 

.1078e-01 

0.30505 

0.29427 

0.29961 

34 

. 1040e-01 

0.29427 

0.28387 

0.28903 

35 

.1003e-01 

0.28387 

0.27384 

0.27881 

36 

•9677e-02 

0.27384 

0.26416 

0.26896 

37 

.93356-02 

0.26416 

0.25483 

0.25946 

38 

.9005e-02 

0.25483 

0.24582 

0.25029 

39 

.8687e-02 

0.24582 

0.23714 

0.24144 

40 

,8380e-02 

0.23714 

0.22876 

0.23291 

41 

.8084e-02 

0.22876 

0.22067 

0.22468 

42 

•7798e-02 

0.22067 

0.21288 

0.21674 

43 

.75236-02 

0.21288 

0.20535 

0.20908 

44 

.7257e-02 

0.20535 

0.19810 

0.20169 

45 

.7000e-02 

0.19810 

0.19110 

0.19456 

46 

•6753e-02 

0.19110 

0.18434 

0.18769 

47 

.651 4e-02 

0.18434 

0.17783 

0.18106 

48 

.6284e-02 

0.17783 

0.17154 

0.17466 

49 

.6062e-02 

0.17154 

0.16548 

0.16849 

50 

.5848e-02 

0.16548 

0.15963 

0.16253 

51 

.5641 e-02 

0.15963 

0.15399 

0.15679 

52 

.5442e-02 

0.15399 

0.14855 

0.15125 

53 

.5250e-02 

0.14855 

0.14330 

0.14590 

54 

.5064e-02 

0.14330 

0.13824 

0.14075 

54 

.4885e-02 

0.1 3824 

0.13335 

0.13577 

56 

.471 2e-02 

0.13335 

0.12864 

0.13097 

57 

.4546e-02 

0.12864 

0.12409 

0.12635 

58 

.4385e-02 

0.12409 

0.1 1971 

0.12188 

59 

.4230e-02 

0.1 1971 

0.11548 

0.11757 

60 

•408le-02 

0.1 1548 

0.11140 

0.11342 

61 

. 3937e-02 

0.11140 

0.10746 

0.10941 

62 

. 3797e-02 

0.1C746 

0.10366 

0.10554 

63 

.3663e-02 

0.10366 

0.10000 

0.10182 
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D.  CALIBRATION  PARTICLES 


The  scatterers  used  for  system  calibration,  theoretical  analysis,  and 
performance  determination  were  polystyrene  latex  spheres  of  known  size 
obtained  from  Duke  Scientific  Corporation  of  Palo  Alto,  California  and  Seragen 
Diagnostics  (formerly  Dow  Chemical  Company)  of  Indianapolis,  Indiana.  These 
particles,  suspended  in  an  admixture  of  water  and  surfactant,  are  listed  in 
Table  3.  Both  companies  report  the  density  as  1  .05  g/mt  and  the  refractive 
index  as  1.59  at  589.26  nm.  Spectron  reports  a  refractive  index  of  1.57  at 
488.0  nm.  All  particle  solutions  from  Duke  Scientific  are  0.5  percent  solid 
by  weight,  except  the  0.412  pm  size  which  is  10  percent  solid  by  weight.  All 

of  the  Dow  solutions  are  10  percent  solid  by  weight.  Samples  are  prepared  by 

withdrawing  solution  from  the  highly  concentrated  solution  with  a  microliter 
syringe  and  diluting  them  with  200  ml  deionized  water  or  methanol.  Particles 
in  this  resulting  solution  are  then  atomized  by  a  TSI  Constant  Output 
Atomizer. 

E.  PARTICLE  GENERATION  SYSTEM 

The  TSI  Constant  Output  Atomizer  (Model  3076)  is  shown  in  Figure  11. 
Compressed  air  is  first  dried  and  filtered  before  being  supplied  to  the 
atomizer  assembly.  In  the  assembly,  the  air  is  expanded  through  a  small  hole 
(0.0135  inches)  to  form  a  high-velocity  jet.  The  solution  is  drawn  by  the  jet 
to  the  atomizer  where  large  droplets  are  removed  by  impaction  on  the  wall  of 

the  atomizer.  The  excess  liquid  returns  to  the  container.  The  aerosol 

travels  to  the  diffusion  drier  (TSI  Model  3062)  shown  in  Figure  12.  The  dried 
aerosols  are  injected  through  a  nozzle  centered  above  the  optical  probe 
volume. 

F.  CALIBRATION  AND  APPLICATION  OF  THE  IPSS 

The  Mie  theory  (References  1,  3 »  and  4)  shows  that  the  scattered  intensity 
ratio  at  angles  of  scatter,  9  and  $,  can  be  found  if  the  refractive  indicies 
of  both  the  particle  and  surrounding  medium,  the  wavelength  and  polarization 
of  the  incident  radiation,  and  the  diameter  of  the  scatterer  are  all  known. 
However,  in  experimental  applications  (soot-sizing  measurements),  the  IPSS 
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TABLE  3-  PARTICLE  SIZES  AVAILABLE  FOR  EXPERIMENTATION. 


MEAN  DIAMETER  (ym)  STANDARD  DEVIATION  ( ym) 


0.087* 

0.0046 

(5.3*) 

0.109* 

0.0027 

(2.5*) 

0.173 

0.0068 

(3.9*) 

0.198 

0.0036 

(1  .8*) 

0.215 

0.0026 

(1.2*) 

0.255* 

0.0022 

(0.9X) 

0.261 

0.0031 

(1.2*) 

0.272 

0.0031 

0.312 

0.0022 

(0.7*) 

0.330 

0.0040 

0.364* 

0.0024 

(0.7*) 

0.412 

0.0058 

(1  .4*) 

0.460 

0.0048 

(1 .0*) 

0.497 

0.0059 

(1 .2*) 

0.500* 

0.0027 

(0.5*) 

0.605 

0.0089 

(1  .5*) 

0.624* 

0.0052 

(0.8*) 

0.760 

0.0046 

0.804* 

0.0048 

(0.6*) 

*  Obtained  from  Duke  Scientific  Corporation 
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igure  11.  TSI  Constant  Output 


requires  the  solution  of  the  inverse  problem— that  is,  the  scattered  intensity 
ratio  is  known  (through  measurement)  and  the  particle  size  must  be 
determined.  As  stated  earlier,  this  correlation  is  provided  by  the  bin  tables 
present  in  the  IPSS  software  package.  However,  to  ensure  that  the  IPSS 
correctly  reports  the  diameter,  the  system  must  first  be  calibrated  by 
injecting  latex  particles  of  known  size  and  adjusting  the  PMT  gains  until  the 
measured  mean  diameter  matches  the  size  of  the  injected  particles. 

Previously,  the  standard  procedure  for  delivering  the  particles  to  the 
probe  volume  for  calibration  was  to  dilute  them  in  deionized  water  and  pass 
the  resulting  suspension  through  an  atomizer  to  form  droplets.  The  droplets, 
then,  are  mixed  with  dry  air,  exposing  the  particles  in  their  natural  size  and 
refractive  index  in  the  probe  volume.  It  was  discovered  that  droplets  of 
water  escaping  from  the  particle-dry ing  system  were  overwhelming  the  signals 
from  the  control  volume.  This  resulted  in  size  measurement  bias  towards  a 
static  distribution;  the  system  could  not  respond  to  different  size 
particles.  Dilution  with  a  more  volatile  non  interfering  liquid,  such  a3 
methanol  allows  the  system  to  respond  effectively  to  different  size 
particles.  The  size  distributions  obtained  were  narrower;  the  repeatability 
of  the  results  was  found  to  be  more  consistent,  and  the  measured  mear;  diameter 
for  each  particle  approached  the  nominal  diameter  specified  by  the 
manufacturer.  Because  alcohols  can  disturb  the  size  of  the  latex  particles 
with  time,  a  new  particle  solution  should  be  used  for  every  calibration  of  the 
instrument.  Figure  13  displays  the  measured  linear  mean  diameter  versus  the 
corresponding  nominal  diameter  of  each  particle  measured. 

During  calibration,  a  particle  solution  was  diluted  to  half  the  original 
concentration  to  check  the  type  of  scattering  occuring  at  the  probe  volume. 
Correspondingly,  the  data  rate  decreased  by  almost  half.  The  above  seems  to 
indicate  that  single  scattering  is  the  most  important  scattering  process 
occuring  at  the  probe  volume.  Figures  14a,  and  14b  show  normalized  size 
distribution  histograms  for  0.412  ym  latex  particle  solutions  of  arbitrary 
concentrations  cl  and  one-half  of  cl,  respectively. 

Also  during  calibration  of  the  40V200  angle  pair,  some  of  the  diameter 
bin  widths  listed  in  the  IPSS  data  output  table  (and  used  in  the  calculation 
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Figure  13.  Measured  Linear  Mean  Diameter  versus  Nominal  Diameter  of 
Latex  Spheres. 


of  the  standard  deviation)  were  found  to  have  negative  bin  widths.  A 

comparison  of  the  theoretical  values  (obtained  from  the  Mle  code  found  in 
Reference  4)  with  the  IPSS  bln  table  data  (intensity  ratio  bin  centers  versus 
diameter  bin  centers  for  the  angle  pair  40°/20°  and  refractive  index  of  latex, 
1.57-iO.O),  shown  in  Figure  15,  revealed  a  disagreement  between  the  two 

plots.  This  disagreement  was  unexpected  since  it  appears  that  the  IPSS  bin 
tables  were  produced  by  an  interpolation  of  intensity  ratios  obtained  from  Mie 
code  output.  The  discrepancy  is  currently  attributed  to  the  assumption  of 
constant  absolute  scattered  intensity  within  the  solid  angle  of  each 
collection  lens.  The  theoretical  plots  shown  in  the  following  discussion  were 
obtained  from  this  assumption  and  are  not  expected  to  be  valid  for  the  solid 
angle  present  in  the  IPSS.  Through  interfacing  the  Mie  code  with  algorithms 
which  will  generate  data  points  on  the  lens  and  integrate  these  points  over 
the  solid  angle  of  the  receiving  lens,  the  theoretical  plot  is  expected  to 
collapse  to  the  plot  provided  by  the  IPSS  bin  tables. 

Although  the  theoretical  plots  do  not  agree  with  IPSS  bin  table  plots,  the 

trends  are  correct.  Hence,  the  presence  of  negative  bin  widths  can  be 

explained.  In  addition,  insight  is  gained  into  the  sensitivity  and 

performance  of  the  instrument.  In  a  closer  examination  of  the  theoretical 
intensity  ratio  curve  at  a  refractive  index  of  1  .57  and  an  angle  pair  of 
40o/20°,  a  local  minimum  and  maximum  between  the  diametric  range  of 
C.53  urn  and  0.60  um  is  found  (Figure  16).  Thus,  in  this  small  region,  the 
intensity  ratio  is  multivalued.  If  such  a  behavior  is  not  considered  when 
developing  and  implementing  an  interpolation  scheme  to  generate  bin  table 
data,  results  such  as  negative  bin  widths  can  be  obtained.  The  bending  of  the 
curve  back  onto  itself  in  the  bin  plot  table,  shown  in  Figure  15,  demonstrates 
the  negative  bin  widths  resulting  from  the  presence  of  a  local  maximum  and 
minimum.  Since  particle  size  cannot  be  resolved  in  these  multivalued  regions, 
calibration  in  or  near  these  areas  must  be  avoided. 

Similar  comparisons  for  the  other  angle  pairs  of  both  the  large  and  small- 
angle  ratioing  technique  can  also  be  made.  As  shown  in  Figures  17  and  1 8 ,  the 
scattered  intensity  ratio  versus  diameter  curves  for  latex  at  the  other  angle 
pairs  of  the  large-angle  ratioing  technique  (50°/20°  and  60°/20°)  do  not 
exhibit  multiplicities  in  the  range  of  measurement  specified  by  Spectron.  It 
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Figure  16.  A  Closer  Examination  of  the  40°/20°  Angle  Pair  for  a  Real 
Refractive  Index  of  1.57  and  the  Diametric  Range  of  0.53 
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Figure  19.  Plot  of  Intensity  Ratio  versus  Diameter  for  an  8°/2.5 
Angle  Pair  and  a  Real  Refractive  Index  of  1.57. 
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Figure  20.  Comparison  of  Theoretical  Intensity  Ratio  versus  Diameter 
with  the  IPSS  Bin  Center  Tables  for  a  10°/5°  Angle  Pair 
and  a  Real  Refractive  Index  of  1.57. 
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i3  expected,  then,  that  no  difficulties  should  be  encountered  when 
implementing  these  angle  pairs.  Unfortunately,  the  same  comments  cannot  be 
made  about  the  small-angle  ratioing  technique  for  latex.  For  such  small-angle 
pairs,  the  presence  of  a  real  refractive  index  is  responsible  for  rapid 
oscillations  in  the  intensity  ratio  versus  diameter  curves  (Figures  19  through 
21).  Again,  the  local  maximums  and  minimums  are  responsible  for  the 
calculations  of  negative  bin  widths.  However,  due  to  the  presence  of  the 
numerous  oscillations,  the  calculation  of  bin  tables  become  difficult,  if  not 
impossible.  The  bin  table  plot  for  the  10°/5°  angle  pair  (shown  in  Figure  20) 
and  the  absence  of  bin  tables  for  the  remaining  angle  pairs  demonstrate  this 
point.  At  its  best,  then,  the  small-angle  ratioing  technique  for  particles 
possessing  a  real  refractive  index  is  severely  limited  to  extremely  small  size 
ranges  outside  the  range  of  oscillations. 

Finally,  for  the  case  of  scatterers  possessing  a  complex  refractive  index 
(Figures  22  through  27),  the  intensity  ratio  versus  diameter  curves  are  well- 
behaved  for  both  small  and  large-angle  ratioing  techniques.  However,  3lnce 
instrument  calibrations  are  carried  out  with  latex  particles  (having  a  real 
refractive  index),  the  "goodness"  of  the  measuring  technique  for  soot  (having 
a  complex  refractive  index)  is  dependent  on  how  well-behaved  the  latex  curves 
are.  Therefore,  the  use  of  the  small-angle  ratioing  technique  for  soot  size 
measurements  does  not  look  promising. 
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Figure  22.  Plot  of  Theoretical  Intensity  Ratio  versus  Diameter  for 
an  8°/2.5°  Angle  Pair  and  a  Complex  Refractive  Index  of 
1.57  -  i0.56. 
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Figure  23 .  Comparison  of  Theoretical  Intensity  Ratio  versus  Diameter 
with  the  IPSS  Bin  Center  Tables  for  a  10°/5°  Angle  Pair 
and  a  Complex  Refractive  Index  of  1.57  -  i0.56. 
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Figure  2  6.  Comparison  of  Theoretical  Intensity  Ratio  versus  Diameter 
with  the  IPSS  Bin  Center  Tables  for  a  50°/20°  Angle  Pair 
and  a  Complex  Refractive  Index  of  1.57  -  0.56. 
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the  IPSS  Bin  Center  Tables  for  a  60°/20°  Angle 
Complex  Refractive  Index  of  1.57  -  i0.56. 


SECTION  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

It  has  been  shown  that  the  scattered  intensity  depends  on  the  magnitude, 
wavelength  and  polarization  of  the  incident  intensity,  the  refractive  indices 
of  the  particle  and  the  medium  outside  the  particle,  the  particle  size  and  its 
geometry,  the  distance  away  from  the  center  of  the  particle,  and  the  angles  of 
scatter.  The  assumptions  made  to  arrive  at  the  Mie  expression  are  as  follows: 

1.  both  e,  u,  and  <j  are  isotropic,  time-independent  electromagnetic 
properties, 

2.  the  electric  and  magnetic  fields  are  spatiall"  and  temporally 
continuous, 

3.  either  the  conductivity,  0,  is  zero  or  is  included  in  the  measurement 
of  the  imaginary  part  of  the  refractive  index, 

4.  the  scattering  is  restricted  to  single,  elastic,  and  independent 
scattering, 

5.  the  particle  is  spherical, 

6.  the  particle  material  is  homogeneous, 

7.  the  magnetic  permeability  of  the  particle  is  equal  to  the  magnetic 
permeability  of  the  outside  medium, 

8.  the  relative  refractive  index — that  is,  the  ratio  of  the  refractive 
index  of  the  particle  to  that  of  the  outside  medium — must  be  known. 

9.  no  free  surface  charge  is  present  on  the  particle, 

10.  the  incident  radiation  is  a  time-harmonic,  plane  wave,  and 

11.  the  scattered  radiation  is  measured  far  from  the  point  of  scatter  so 
that  the  longitudinal  electric  field  is  negligible. 

The  IPSS  design  complies  with  the  last  two  limitations  by  the  placement  of 
the  detectors  in  the  far  field  and  the  implementation  of  the  laser  as  the 
incident  source.  The  remaining  limitations  are  imposed  on  the  system  by  the 
shape  and  composition  of  the  particle  and  are,  therefore,  possible  sources  of 
uncertainty  in  measurement.  The  effects  of  particles  with  varying  refractive 
indices  and  aspherical  shape  on  the  accuracy  of  an  instrument  using  a 
variation  of  the  small-angle  ratioing  technique  have  been  studied  by  Reference 


8.  Although  the  conclusions  drawn  in  that  study  may  not  be  valid  for  the 
large-angle  pairs  of  the  IPSS,  the  same  experimental  procedures  can  be  used  to 
test  the  IPSS's  response  to  these  effects.  The  effect  of  charged  particles  on 
the  system  is  also  a  topic  for  further  research  since  these  particles  are 
present  during  combustion.  Reference  9  has  presented  the  groundwork  for  the 
theoretical  study  of  this  problem. 

In  Section  III,  it  was  proposed  that  the  theoretical  plot  would  collapse 
to  the  IPSS  bin  table  plots  once  integration  over  the  solid  angle  was 
performed.  Preliminary  results  have  shown  that  this  is  not  the  case.  In 
fact,  the  integrated  intensity  values  do  not  differ  significantly  from  the 
theoretical  point  values  plotted  in  Figures  15  and  IT  through  27.  Neither  a 
correlation  nor  an  explanation  has  been  found  which  explains  this 
disagreement.  In  this  light — assuming  no  oversight  has  been  committed  in 
obtaining  these  results — the  IPSS  tables  for  the  large-angle  technique  should 
be  replaced.  However,  the  ability  to  generate  new  tables  also  provides  the 
ability  to  make  the  IPSS  more  versatile  since  parameters  such  as  refractive 
index  and  angle  pair  can  now  be  varied. 
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APPENDIX  A 


MAXWELL’S  EQUATIONS 


The  basic  laws  of  electromagnetic  radiation  are  known  as  Maxwell's 
equations  and  are  listed  below: 


(A-1  ) 
(A-2) 
( A— 3 ) 
(A-4) 

where  D-  electric  flux  density  or  displacement  flux  density  [C/m2], 

E  »  electric  field  intensity  [volts/m], 

B  »  magnetic  flux  density  [Webers/m2], 

H  =  magnetic  field  intensity  [A/m], 
p  =  charge  density  [C/m^], 

-4  p 

J  =■  conduction  and  convection  current  density  [A/m  ]. 


7  •  D 
7  •  B 
7  x  E 


3B 

at. 


V  X  H  -  j  /§£ 


Gauss's  Law  for  Electric  Fields 
Gauss's  Law  for  Magnetic  Fields 
Faraday's  Law 
Ampere's  Law 


These  terms  are  also  related  by  the  following: 

(A-5) 
( A-6 ) 
( A— 7 ) 

where  e  =  electric  permittivity  [farad/m], 
p  =  magnetic  permeability  [Henry /m], 
o  «  conductivity  [mhos/m]. 


D  =  eE 

-*  -¥ 

B  =  pH 

-♦  4 

J  =  oE 


Although  the  point  form  of  Equations  (A-1)  to  (A-1*)  may  appear  confusing 
and  meaningless  to  those  who  have  not  studied  electromagnetic  field  theory,  a 
physical  representation  can  be  drawn  from  the  integral  form  of  the  equations. 
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A.  GAUSS'S  LAW  FOR  ELECTRIC  FIELDS 


By  integrating  7  •  D  *  p  over  any  volume, 


J  ( 7  •  D)dV  =  /v pdV  =  Q 


where  q  »  total  free  charge  in  the  volume,  and  applying  the  divergence 
theorem, 


[  (7  •  F)dV  •  F  •  dS  , 

1 V  ■'3 


the  integral  form  can  be  obtained 


P. 


This  equation  states  that  the  amount  of  charge  enclosed  by  a  surface  can  be 
related  to  the  lines  of  flux  entering  (for  negative  charges)  or  leaving  (for 
positive  charges)  a  surface  (See  Figure  A-1). 

B.  GAUSS'S  LAW  FOR  MAGNETIC  FIELDS 

In  a  similar  manner,  7  •  B  =  0,  integrated  over  a  volume  becomes 


=  0 


This  equation  states  that  the  net  magnetic  flux  through  a  closed  surface  is 
always  zero.  In  other  words,  the  magnetic  flux  lines  close  upon  themselves 
(See  Figure  A-2).  The  magnetic  counterpart  of  an  isolated  electric  charge — 


Figure  A-l.  A  Physical  Representation  of  Gauss's 
Law  for  Electric  Fields 


namely,  a  single  magnetic  pole — does  not  exist.  For  every  north  pole,  a  south 
pole  exists. 

C.  FARADAY'S  LAW  OF  INDUCTANCE 

If  the  point  form  of  Faraday's  Law  is  integrated  over  a  surface, 


j  (7  x  E)  •  dS  =  -  j  *  dS 
J  3  •'sot 


Stokes  Theorem  can  be  applied  to  the  left  hand  side  to  obtain 


J  E 


•  dl 


f  3B 

"  J  at 


dS. 


It  can  be  seen  that  a  time  varying  magnetic  field  (H  «  B)  induces  an 
electromotive  force 


emf  =  <£  E  •  d& 


which  establishes  a  current  in  a  closed  circuit.  An  example  of  such  a  case  is 
shown  in  Figure  A-3.  Here  the  movement  of  the  magnet  toward  the  closed 
circuit  causes  an  increase  in  the  magnetic  flux  through  any  arbitrary  surface, 
S,  enclosed  by  the  coil,  and  as  a  result,  the  potential  drop  across  the 
resistor  increases  in  such  a  way  as  to  oppose  the  increasing  magnetic  field. 
Note  that  although  the  previous  example  employed  a  physical  closed  conducting 
path,  the  path  could  be  an  imaginary  line  in  space. 

D.  AMPERE'S  LAW 

The  integral  form  of  Ampere's  Law  can  be  written  by  again  employing 
Stokes'  Theorem  to  obtain 
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j  H  •  d£  -  Jgj  •  dS  +  /  |~  •  dS. 


By  dimensional  analysis,  it  can  be  seen  that  the  right  hand  side  of  the  above 
equation  are  current  terms  where  the  first  term  is  the  conduction  or 
convection  current,  and  the  second  term  is  the  displacement  current.  For 
example,  consider  the  circuit  shown  below  in  Figure  A-4.  If  a  time-varying 
magnetic  field  pierces  the  plane,  S,  of  the  circuit,  by  Faraday's  Law,  an 
electric  potential  is  induced  in  the  conductor  causing  current  to  flow.  The 
integration  of  the  magnetic  field  tangential  to  the  path,  da,  yields  the 
current  exiting  the  arbitrary  surface  enclosed  by  the  path.  At  (a)  the 
current  leaving  the  surface  is  the  conduction  current 


and  the  displacement  current. 


r  3D 
JS  3t 


dS 


is  zero.  If  the  surface  is  made  to  cut  between  the  plates  of  the  capacitor  at 
(b),  Ic  =  0  because  charge  is  not  being  transmitted  from  one  plate  to 
another.  However,  the  line  integral  still  follows  the  same  path;  therefore, 
the  total  current  calculated  by  each  method  must  be  equal.  Thus,  a 
displacement  current  exists  between  the  plates  of  the  capacitor  causing  the 
circuit  to  be  complete. 

E.  BOUNDARY  CONDITIONS 


The  previous  equations  will  give  expressions  of  the  electromagnetic  field 
as  it  propagates  through  a  medium.  However,  if  another  medium  is  encountered, 
boundary  conditions  must  be  known  at  the  interface. 
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From  the  integral  form  of  Faraday's  Law  and  Ampere's  Law  respectively,  it 

can  be  shown  that  the  electric  field  and  magnetic  field  components  tangential 

to  the  surface  of  the  interface  are  continuous  provided  that  J,  and  — 

ot  3t 

remain  finite  as  the  surface  area  approaches  zero.  Thus, 


Et,1  =  Et,2 


H 


t,1 


\,2‘ 


(A-8a,b) 


From  Gauss's  Law,  the  normal  components  of  the  electric  flux  density,  and 
therefore,  the  electric  field,  are  in  general  discontinuous  across  the 
boundary.  Thus, 


Jn,  1 


-  D, 


n  ,2 


or 


e,E  .  ■ 
1  n  f  1 


£2En,2 


=  P. 


(A-8c) 


where  pg  is  the  free  surface  charge  density  existing  on  the  boundary. 
Applying  the  same  law  to  the  normal  magnetic  field  components,  it  is  seen  that 


B 

n, 


1 


or  m,H  , 
1  n,  1 


2  n ,  2 


(A-8d) 


If  the  nature  of  the  medium  (i.e.,  conductor  or  insulator)  is  known,  the 
boundary  conditions  can  be  made  more  specific.  Consider  first  the  case  of  the 
perfect  conductor.  Because  no  electromagnetic  fields  exist  in  the  interior  of 
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a  perfect  conductor  (Figure  A-5),  the  boundary  conditions  become 


(A-9a,b,c,d) 


where  K  is  the  current  per  unit  width  flowing  on  the  surface  of  the  conductor 
perpendicular  to  the  tangential  magnetic  field  at  the  boundary  of  the 
surface.  Note  for  a  perfect  conductor,  that  although  the  current  flowing  on 
the  surface  is  finite,  the  current  density  J  is  infinite.  Therefore,  Equation 
(A-8b)  is  no  longer  valid. 

For  a  perfect  insulator  (dielectric),  no  free  surface  charge  is  present  on 
the  surface  unless  the  charge  was  specifically  placed  there.  In  the  case  of 
no  free  surface  charge,  the  boundary  conditions  become 


n,  1 

=  °n,2 

:t,1 

=  Et ,  2 

[t,1 

“  Ht  ,2 

'n,1 

*  Bn ,  2 

(A-1 Oa , b , c , d ) 
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